QUASI-INVARIANT MEASURES, ESCAPE RATES 
AND THE EFFECT OF THE HOLE 
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Abstract. Let T be a piecewise expanding interval map and Th be an ab- 
stract perturbation of T into an interval map with a hole. Given a number I, 
< £ < 1, we compute an upper-bound on the size of a hole needed for the 
existence of an absolutely continuous conditionally invariant measure (accim) 
with escape rate not greater than — ln(l — i). The two main ingredients of 
our approach are Ulam's method and an abstract perturbation result of Keller 
and Liverani. 



1. Introduction 

Open dynamical systems have recently been a very active topic of research in er- 
godic theory and dynamical systems. Such dynamical systems are used in studying 
nonequilibrium statistical mechanics [19] and metastable chaos [22|. 

A dynamical system is called open if there is a subset in the phase space such 
that whenever an orbit lands in it, the dynamics of this obit is terminated; i.e, 
the orbit dies or disappears. The subset through which orbits escape is called 
a hole, denoted H . The escape rate through H can be measured if the system 
admits an absolutely continuous conditionally invariant measure (accim). The first 
result in this direction is due to Pianigiani and Yorke [15] • The survey article [S] 
contains a considerable list of references on the existence of accim and its relation 
to other measures. One of the most intuitive existence results is found in Section 
7 of |16j . It is mainly concerned with systems having small holes and its idea is 
based on the perturbation result of [lOj . It roughly says that if a mixing interval 
map is perturbed by introducing a 'sufRciently small' hole, then the resulting open 
dynamical system admits an accim. Our main goal in this paper is to show how 
the condition 'sufRciently small' can be computationally verified in some of these 
results, in particular, results from Section 7 of [16 . 

More precisely, for a given Lasota- Yorke map T, we use Ulam's method on the 
closed dynamical system T to give a computable size of the hole H for which the 
open dynamical system T}j must admit an accim. 

Historically, Ulam approximations have been used to provide rigorous estimates 
of invariant densities of closed systems (see [15] and references cited there) or to 
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approximate other dynamical invariants (see 0). The use of Ulam's method in 
the study of open systems is natural. In [1], Ulam approximations were used to 
rigorously estimate the escape rate for certain open systems. In that computation 
the Ulam matrix was derived from the Perron-Frobenius operator associated to 
the open system (i.e., a sub-stochastic matrix), whereas here, we approximate the 
closed system. The method in [Tj also demanded existence of an accim as a basic 
assumption. 

We remark that, as a consequence of the spectral methods discussed here, upper 
bounds on the escape rate can be obtained from analysis of the closed system. 
However this does not generally replace the computation in [1] where the Ulam 
approximation of the open system yields both upper and lower bounds on the 
escape rate (but under the additional assumption of existence of an accim). Hence, 
there is potential to apply a two step method - the current algorithm would be used 
to guarantee an accim and to provide rough (upper) bounds on the escape rate, 
followed by the method of , once the size of the hole is fixed, to more accurately 
estimate the latter. 

Our paper is organized as follows. In Section 2 we present a version of Keller 
and Liverani's abstract perturbation theorem. The constants which are involved 
in this theorem are essential in all our computations and thus, we need to state 
all the details of this theorem explicitly. Section 3 contains a precise setting of 
the problem. Section 4 contains technical lemmas, mostly well-known and stated 
without proof. Section 5 presents the algorithm (Algorithm 15. ip whose outputs 
(Scorn, £ com) are the parameters used to solve our problem. Theorem l5.31 takes these 
parameters and computes a maximum hole size leading to an accim. In Section 6 
we provide, in detail, rigorous computations of the size of a hole for two examples 
as well as a discussion of computational overhead and some techniques for reducing 
computation time. In Section 7 we discuss how our methods can be implemented 
in a smooth setting where there are interesting results concerning the effect of the 
position of a hole on the escape rate. This is in connection with the recent results 
of [1] and [Tl] . The examples of Section [S] and the discussion of Section [7] highlight 
the new results that Algorithm 15. II brings to open dynamical systems. 

2. The First Keller-Liverani Perturbation Result 

Let (/, 05, A) be the measure space where I — [0, 1], 03 is the Borcl a-algcbra and 
A is Lebesgue measure. Let = L^{I, S, A). For / G i^, we define 



We denote by BV the space of functions of bounded variation on I equipped with 
the norm || ■ \\bv ^ Vi-) || • ||i [7]. Let : BVil) BV{I) be two bounded 
linear operators, i — 1,2. We assume that: For f £ 



Vf = inf{var/ : / = / a.e.}. 



/ 



where 



i-i 




(2.1) 



ll^./lli < ll/lli 
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and 3 a G (0, 1), A > and B >0 such that 

(2.2) ||i'r/llBV<^«"||/||w + i?||/||i VneN V/e ^ = l,2. 

Further, we introduce the mixed operator norm: 

|||P,;|||= sup llPJIli. 

II/IIbv<i 

For any bounded hnear operator P : BV — > BV with spectrum <j{P), consider the 
set 

VsAP) = {z e C : |z| < r or dist(z, cr(P)) < 5}. 

Since the complement of Vs^r{P) belongs to the resolvent of P, it follows that ([7] 
Lemma 11, VII.6.10) 

z 

Remark 2.1. a in (|2.2p is an upper bound on the essential spectral radius of P^ [2]. 

Theorem 2.2. [TUIITQ Consider two operators Pi : BV{I) BV{I) which satisfy 
|gj]) and (KM>- For r G (a, 1), let 

ni = \-, ^1 

mr/a 

C^r-''^; D^A{A + B + 2) 

An%BDCHsAPi). 

"2 = I -. ■ 

in r / a 

If 

„ni+n2 , , 

then for each z G C \ VsAPi)j ™s have 

\\{z ~ P2)-'f\\BV < ^^r-^ WfWBV + ^ll/lli. 

ln(r/ck) 
" ln(l/a)' 
8[2A{A + B) + {1- r)-^]{A + Bfr-^^ + 1 



Set 



and 



1 -r 

b=2[{A{A + Bf{D + B) + B){l-r)-\-'"'^ +P]. 



4P {HsAPi){D + B) + 2A{A + B) + {1 - r)-i) 



If 

(2.3) 

IIIP1-P2III <min{ei(Pi,r,(5), 

=^eo(Pi,r, S) 

then for each z G C \ r{Pi), we have 
(2.4) 

|||(z-P2)-i-(z-Pi)-i|||<|||Pi~P2||r(a||(z-Pi)-i||BV + 6||(z-Pi)-ily) 
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Corollary 2.3. [IHl [ig //|||Pi - P2III < Si{Pi,r,S) then (7(^2) C Vs^riPi)- In 
addition, i/ |||Pi — P2III < ^oiPiTf, , then in each connected component ofVs^riPi) 
that does not contain both cr(Pi) and cr{P2) have the same multiplicity; i.e., the 
associated spectral projections have the same rank. 

3. Expanding Interval Maps, Perturbations and Holes 

Let r be a non-singular interval map and denote by P the Perron-Frobenius 
operator associated with T [5]. Typically and this will be the case for our examples, 
T will be a Lasota-Yorke mapQ- We assume: 

(Al) 3 ao e (0, 1), and Pq > such that V/ G BV{I) 

VPf <aoVf + Bo\\f\\i. 

Remark 3.1. Condition (Al) implies that p = 1 is an eigenvalue of P. In particular, 
T admits an absolutely continuous invariant measure [21 [3]. Moreover, for any 
r G (ao, 1) there exists a So > 0, 60 depends on r, such that for any Sq G (0, 60] and 
any eigenvalue pi of P, with \pi\ > r, we have: 

(1) P(p„(5o)nP(0,r) = 0; 

(2) B{p,,So)nB{p,,5o)=9,i^j. 

The inequality of assumption (Al) is known as a Lasota-Yorke inequalitj0. For a 
given Lasota-Yorke map with (3 > 2, the constant ao 2//3 in inequality (Al) and 
Po may be found in terms of bounds on the second derivative of the map T and 
the minimum of x^-j-i — Xi. 

3.1. Ulam's approximation of P. Let rj he a finite partition of / into intervals. 
Let mesh(77) be the mesh size of 77; i.e, the maximum length of an interval in 77, and 
let be the finite a-algebra associated with rj. For / G P^, let 

n,/ = E(/|'B,), 

where E(-|5B,,) denotes the conditional expectation with respect to S,,. Specifically, 
if x G G 77 



(n,/)(x) = 




Define 

p,, = n^oPon^. 

Pri is called Ulam's approximation of P. Using the basis {j^~jXin} can 
be represented by a (row) stochastic matrix acting on vectors from M'*^'') by right 
multiplication: x — > xPrj. The entries of Ulam's matrix are given by: 

^ A(/,nr-V,) 



^ A map T acting from an interval [a, b] to itself is a Lasota- Yorke map if it is piecewise 
with respect to a finite partition a = xq < xi < ■ ■ ■ < x„ = b, has well-defined left and right limits 
of derivatives up to second order at each Xi and is expanding: /3 ;= inf^^^^ |T'(a;)| > 1. 

^In fact, (Al) is slightly stronger than the original Lasota-Yorke inequality. In particular, 
when T is a general piecewise expanding map the Lasota-Yorke inequality is given by VPf < 
2/3~^y/-|-_Boll/l|i- El for details and for generalizations of the original result of I14| . In certain 
situations, in particular, when T is piecewise expanding and piecewise onto or when infa, |T'(a;)| > 
2, the original Lasota-Yorke inequality reduces to (Al). In principle, when dealing with Ulam's 
scheme for Lasota-Yorke maps, as we do in this paper, (Al) cannot be relaxed. See [17] for details. 
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Since P,, can be represented as stochastic matrix, it has a dominant eigenvalue ~ 
1 [15 • Any associated left eigenvectors represent invariant functions G BV(V) 
for the operator P,,. 

3.2. Interval maps virith holes. Let iJ C / be an open interval. Denote by 

dcf 

Th := T|Xo, where Xo = I\H. We cah Th an interval map with a hole; H being 
the hole. Its Perron- Frobenius operator, which we denote by Ph, is defined as 
follows: for f E and n > 1 

where X„_i = n"Jg^T~*Xo, the set of points whose orbits do not meet the hole H 
in the first n — 1 steps. 

Definition 3.2. A probability measure /i on [0,1], dji = is said to be 

an absolutely continuous conditionally invariant measure (accim) if there exists 
< e// < 1 such that PhIh = ^hIh- In this case — Ine/f is the escape rate 
associated to fi. 

4. Some Technical Lemmas 

We state (for the most part, without proof) some well-known technical results 
to be used in our computations later in the paper. 

4.1. Lasota-Yorke inequalities and estimates on the difference of opera- 
tors in the mixed norm. Lasota-Yorke inequalities for P and are standard 
results from the literature (see [20] for an original source) . The inequality for Ph is 
not so well-known, but it is straightforward and we derive it below for completeness. 

Lemma 4.1. The operators P and Pri satisfy a common Lasota-Yorke inequality 
as follows: Vn G N V/ £ BV{I) 

\\PV\\BV<a-\\f\\BV + B\\f\\i 

\\P:;f\\BV<alS\\f\\BV + B\\fh 
with B = 1 + i^ao ' ^Independent of rj. For T^, the map with a hole, and under the 
stronger assumption ao < 1/^ we have: Vn £ N V/ G BV{I) 

\\P'Af\\BV<a^f\\Bv + B\\f\\,, 

where a ~ 3ao < 1 and B ~ ^""^gf" with constants independent of H . 

Proof. For the operator Ph, by assumption (Al), for / G BV{I) we have, 

VPnf - VPifxxo) < aoVifxxo) + So||/||i 

<ao(V/ + 2sup/(a;))+ Poll/Ill 
(4.1) ^e-f 

<ao(V/ + 2y/ + 2||/||i) + Poll/Ill 
= a^/ + (2ao + So)||/||i. 



The assumption that ao < 1/3 can be relaxed simply to ao < 1- This relaxation can still 
produce a common Lasota-Yorke inequality for P, and Ph with constants A and B independent 
oi H. A common Lasota-Yorke inequality for P and Ph with ao < 1 can be found in section 
7 of |16) . Hence the results of this paper are still valid for 1/3 < ao < 1 using the appropriate 
Lasota-Yorke inequaility. For the purpose of the examples which we want to present in Section |6] 
it is more sensible to use the constants produced by Lemma 14. II 



6 



Quasi-Invariant Measures, Escape Rates And The Effect Of The Hole 



Therefore, 

n 

VP^J < a^Vf + a'=-i)(2ao + 
fc=i 

and consequently, for all n > 1 

\\PH.f\\BV<a^f\\BV + B\\f\U. 



□ 



Remark 4.2. Since ag < a and _B < _B in the above lemma, we can obtain a 
common Lasota-Yorke inequality, independent of H and rj using coefficients a and 
B: 

\\P:f\\BV<a^f\\BV + B\\fh 
where represents any of the three operators under discussion. For ease of exposi- 
tion, we will use this common inequality in what follows. (However, see Example l6.2l 
in Section [6] for a discussion indicating how using the strongest possible inequality 
can significantly reduce computational overhead.) 

Lemma 4.3. Let T = maxjao + 1, Bq} and e — mesh^rj). 

(1) <re. 

(2) If\{H) < Te then - Ph||| < 2Te. 

Proof. The first statement is standard. For the proof of the second statement, let 
/ G BV{I) and observe that 

m - PH)f\\i < m - n/iii + \\{p- 

< ePll/llw + KH)\\f\\BV < 2er||/||w. 

□ 

4.2. Computer-assisted estimates on the spectrum of P. All the constants 
arising in Theorem 12.21 are (in principle) computable for the finite-dimensional 
operator P^. Thus, as proposed in [TS], we are going to apply Theorem 12.21 with 
Pri as Pi and P as the perturbation P2. This entails some a priori estimates. 

Lemma 4.4. Given P, S > and r e (a, 1), there exists S2 > such that for each 
rj with < mesh{r]) < 82, we have 

(4.2) mesh{rj)<{2T)-^eo{P^,r,5), 
and 

(4.3) |||P, -P||| < i£o(P^,r,<5). 

Proof. See Lemma 4.2 of [H]. □ 

The computation of a lower bound on e(){Pri, 6, r) follows the argument in Lemma 
3.10 of Ij the only difference arising from the fact that here, P is associated to T 
and not Th as in [T|. The key idea is to estimate the BV-norm of the resolvent of 
Pfj (difficult to compute) by the || • ||i— norm. Hence, following Lemma 3.10 of [1] 
we define 

i?|.,(P„) = sup{(-^ + l)||(z-P^)-i«|li + ^^ + - : ll^lli = l,z e C\n>(P^)}, 
r ~ ao r — ao r 



^In fact, more precisely, we make our computations on the matrix representation of P,, acting 
on the basis {^^^7— yX/, } ^ Section 3. 
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8BiH*{P^)B + (1 - r) 



AB[H*JP^){D + B) + 2(1 + B) + (1 - r)-i 
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Lemma 4.5. 

(1) eQ{Pri,6,r) is uniformly bounded below; 

(2) e*„iPr„6,r) <eo{PrnS,r); 

(3) mesh{r]) < (2r)^^eQ(P^, 5, r) implies mesh{ri) satisfies J^.iip . 

Proof. Follow the proof of Lemma 3.12 [1] verbatum. □ 

5. Main result 

Now we have our tools ready to use the computer and rigorously solve the fol- 
lowing problem: Given a map T satisfying (Al) and given a number i > compute 
a number e > such that if A(iJ) < Te then the map Th has an accim with escape 
rate — In e/f < — hi(l — £). 

The critical step is to obtain control on the separation of the point spectrum of P 
outside the essential spectral radius a. Naturally, from a computational viewpoint 
we can only really do this for after which we use Theorem 12.21 to transfer the 
picture to the spectrum of P. 

More precisely, the following algorithm will, given the number £ with < £ < 
1 — a, compute a number S — Scom with < S < £ and e — Scom > such that with 
r = I — £, and any rj with mesh(r/) < e 

(1) mesh(77) < (2r)-i£o(P^,r,(5); 

(2) _B(1, S) n B{pi,6) — 0, whenever pi is an eigenvalue of P,,, pi ^ P(l, S) and 

\pi\ > r. 

Thus we obtain the required spectral separation (near the eigenvalue 1) for P^ as 
well as the conditions necessary to apply Theorem [ 



Algorithm 5.1. T and £ given as above, then 

(1) Setr^l~£. 

(2) Pick S ^ ^ < £, k eN. 

(3) Feed in a partition of I into intervals. Call it rj. 

(4) Compute e the mesh size ofrj. 

(5) Find Pjj = (Pi^j,^) where 

A(/,nr-ij,) 



(6) Compute the following: Hl^{Pr,), = Ti^l, C = r-"i, D = i + B, 

"2 - I n:77^ 1 . 7 - i„(i/a) > -tf - i_„ , i - maxji + ao, -Ool- 

(7) Check if e <(2T)-^el(Pr^,5,r). 

If (7) is not satisfied, feed in a new r] with a smaller mesh size and repeat 
(3)-(7); otherwise, continue. 



(8) List the eigenvalues of Pri whose modulus is bigger than r: pri^i, i — 1, . . . , d. 
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(9) Define: 

CL = { all the eigenvalues from the list which are in B(l, 6)}. 

(10) Check that if p,,,, ^ CL, then B{p,^.„ 5) Ci B{1, S) = 0. 

(11) // (10) is satisfied, report Scom ■= S and £com. '■— e; otherwise, multiply k by 
2 and repeat steps (2)-(ll) starting with the last r] that satisfied (7). 

Proposition 5.2. Alaorithm \5.1\ stovs after finitely many steps. 

Proof. By Lemma for each 6 > and rG(a,l)3e = mesh(7y) > such that 

e<i2Tr'£*iP^,r,6). 

Therefore, the internal loop of algorithm 1 5. II (2)-(7) stops after finitely many steps. 
To prove that the outer loop of stops after finitely many steps, observe that there 
exist a K G N, K < +oo, such that 6 = ^ < inm{£,So}, r = 1 — i and rj with 
£ — mesh(?7) > such that 

e < min{(2r)-ieo(P^, r, S), i2Ty^soiP, r, S)}. 
This implies cr{P,-i) C Vs^r{P) C Vj^ ri-^)- Thus, any P,j eigenvalue which is not in 
CL is contained in B{0, r) or it is at distance of at least S from B{1, S). By Remark 
Em (11) of Algorithm ED is satisfied for this K. □ 

5.1. A computer assisted bound on a hole size ensuring the existence of 

ACCIM. Given the output £com and Scom from Algorithm lS.ll it is now straightfor- 
ward to prove the existence of an accim for Th. As a byproduct of the computation, 
the spectral information obtained from the algorithm shows that the associated es- 
cape rate is at most — ln(l — £). 

Theorem 5.3. Let Th be a perturbation of T into an interval map with a hole. If 
\{H) < re com then: 

(1) Ph has dominant eigenvalue en > whose associated eigenfunction > 
is the density of a Th -accim; 

(2) 1 ~ en < Scom.; 

(3) l-eH<{l+'-^)X{H). 

Proof. Let X{H) < Fecom and set mesh(?7) — Scom- Then by (2) of Lemma [4.31 we 
have 

1 1 ^„ - ^'ff 1 1 1 < 2reeom <eo{P,„l~i,S 

com / • 

Using Corollarv l2 . 31 with P,f = Pi and Ph — P2 we obtain that o-{Ph) C Vs^^^^riPi-;)- 
Now, from Algorithm 15.11 recall that B{l,Scom) H B{0,r) = and if \pi\ > r, pi 
is an eigenvalue of F,,, is not in CL, then B{l,Sconi) H B{pi,5com) ~ 0- Then 
CoroUarv 1 2 . 31 implies that the spectral projections of and Ph on B{1, (5com) have 
the same rank. Hence, Ph must have at least one isolated eigenvalue in B{1, (5com)- 
Let ch denote the spectral radius of Ph. Since Ph is a positive linear operator, 
ch e (j{Ph). Moreover, Ph has isolated eigenvalues in B{l,5com)- Thus, en is an 
eigenvalue of Ph and it must be in B{1, ^com)- This ends the proof of the first two 
statements of the theorem. To prove (3) of the theorem, we first find a uniforn0 
upper bound on the BV-novm. of By Lemma |4. II we have 

V{eHf*H) = VPnfh < o^Vf*H + {2ao + Bo)\\f*H\\i- 



By uniform we mean here an upper bound which is independent of H. Hence it holds for all 
Th with X{H) < recom. 



W. Bahsoun And C. Bosc 



9 



Therefore, 

2ao + Bn 2ao + Bp 

V Jh ^ ~, 'o 1 

eu — a 1 — I — a 

and hence we obtain 

\\fH\\BV<l+ 

Using the fact that P, the Pcrron-Frobcnius operator associated with T, preserves 
integrals, we obtain 



(5.1) 



I C PrndX-eH r fhdXl 
Jq Jo 

I / P/^dA - / PHflidXl 
Jo Jo 



<\\\p-PH\\\-\\rH\\Bv<{i + ^^^^)xm 

I — I — a 

□ 

6. EXAMPLES 



In this section we implement Algorithm 15.11 and Theorem 15.31 of the previous 
section on two sample computations. Our aim is to show the feasibility of the 
computation, while at the same time, to discuss some analytic techniques that can 
be used to reduce the weight of computations for some of the larger matrices P,, that 
may arise during application of Algorithm 15. II Large matrices should be expected 
when a is close to 1 or alternatively, when the escape rate tolerance I is small. We 
will take advantage of the second mechanism; in both examples we use the same 
map T: 

9^ for < a; < ^ 



T{x) 



lOx -i ioT ^ < X < ^ 



10 ^ — 10 

where i = 1,2, ... ,9. However, in the first example £ — 1/25 and in the second 
example I = 1/40; i.e., in second example we will be looking for the size of a hole 
which guarantees the smaller escape rate. We remark none of our computations 
are particulary time consuming except for the computation of an upper bound on 
i?|^(P,,)Q. We now turn to the computations. 
The Lasota-Yorke inequality for P is given by: 

^P/< l/9F/ + 2/9||/||i. 

Therefore P satisfies (Al) with ao = 1/9 and Bp = 2/9 and consequetly, F = 10/9, 
a = 1/3, B ^ 5/3 and D = A{A + B + 2) = 14/3. 

Example 6.1. Given £ — 1/25, using Algorithm 15.11 we show that if X{H) G 
{0, fx 10-% Th has an accim with escape rate — Ine/f < — ln(24/25). The values 
of the variables involved in the computation are summarized in Table 1. 



''In particular, creating an Ulam matrix of size 5000 X 5000, or even much bigger, is not really 
time demanding. Once a computer code is developed for this purpose, which only requires the 
formula of the map and the number of bins of the Ulam partition as an input, it will excute the 
nonzero entries in few minutes if not less. 

^In our computations we found a rigorous upper bound on (P,,) for an Ulam matrix of 
size 5000 X 5000. This computation took few hours using MATLAB on a desktop computer. 

^AU these variables depend on r and 5. Hg ^ and n2 also depend on e = mesh(j7) . 
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We present here the method which we have foUowed to rigorously compute an 
upper bound on H*g^{P^), for mesh(7?) = 2 x 10"^. Using MATLAB we found 
the dominant eigenvalue 1 of P,, is simple and that there are no other peripheral 
eigenvalues. Moreover, the modulus of any non-peripheral eigenvalues is smaller 
than ao = 1/9. Therefore we have the following estimate (see [9]) 

(6.1) ||(z-P,)-i||i<<5-i||ni||i + ||i?(z)||i, 

where ||ni|| is the projection associated with the eigenvalue 1 of the operator P^, 
and R{z) is the resolvent of the operator P,,(l — IIi). Since |z|>r=l — ^>q;o, 
R{z) can be represented by a convergent Neumann series. Indeed, we have 

(^.(1-ni))", 



(6.2) 




= 7.444310493. 



The computation of the estimate in (16. 2p is the most time consuming step in the 
algorithnij. Using the definition of and inequality ()6.ip . we obtain that 



Hl^ < 45.46070939. 



r 


24/25 


5 


1/26 


e 


2 X 10"* 


Hg^r 


45.46070939 


m 


1 


C 


25/24 
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(2r)-^£5 


0.0002319492040 


Loop I 


Pass 


Loop II 


Pass 


Output I 


Ecom = 2 X 10-^<5com = 1/26 


Output II 


X{H) e (0, ^ X 10"*] =^ Th admits an accim fi 
with escape rate — In en < — In(24/25) 



Table 1. The output of Algorithm O for / = 1/25 



Example 6.2. Given £ — 1/40, using Algorithm 15.11 we show that if X{H) S 
(0,^ X 10"^], Th has an accim with escape rate — Ine// < — ln(39/40). The values 
of the variables involved in the computations are summarized in Table 2. 



Precisely, the work is in the computation of the powers (P7j(l — ni))", n = 1, . . . , 6. Once 
these powers are known the computation of the norm is very fast. However, we will see in the 
next example how we can benefit from these numbers and avoid time consuming computations 
when dealing with a higher order Ulam approximation in the case of a smaller hole. 
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Here, we explain how we have obtained some of the values which appear in Table 
2. In particular, we will explain how we have avoided time demanding computation 
of Hg ^{Pri') and rigorously estimated 

HsAPr,') < 1036.693385, 

where mesh(77') < mesh(77) = 2 x 10"''. In the first pass through the Algorithm lS-H 
we start with e — mesh(7y) = 2 x 10~^. This is the same e which closed Algorithm 
IS.ll in Example 1 6. II So the numbers ||(Pr,(l — ni))"||i, n = 1, . . . , 6, can be obtained 
from the computation of Example l6.1l as they do not depend on r and S. Thus, the 
rigorous estimate 

HlriP^) < 63.73181657 

easily follows. However, the inner loop of Algorithm 15. II will fail because 

2 X 10""* = e = iaesh{ri) > (2r)-ie* = 0.0001763820641. 

Next, Algorithm 15.11 asks us to feed another Ulam partition 77' with mesh(77') < 
2 X 10~^ and to repeat the inner loop of Algorithm l5.1l Here, we have used a 3-step 
trick to avoid time demanding estimate of the new value i7|^(P^/): 

(1) Let us suppose for a moment that we are only concerned with a rigorous 
approximation of the spectrum of P, the operator associated with T. Then 
P, Pjj and Pj^i satisfy a common Lasota-Yorke inequality which does not 
involve a, but rather ao (see Remark l4.2p . Now, re-checking the computa- 
tions which were obtained in the first run of Algorithm 15.11 and this time 
with a = ao — 1/9 and the modfications of B to B = 1 + and D to 
D ^3 + B. Then the value of {2r)-^e^ changes to 

(2r)~i)eo = 0.0002425063815 > 2 x 10"* mesh(?7). 

Consequently, for any r/ with mesh(?7') < mesh(?7), we have 

< rmesh(r7) + rmcsh(?7') 

< 2rmesh(77) < e*{Pr„r,S). 

Therefore, we can use part one of Theorem 12 . 21 with Pi = Prj and P2 = Pry ■ 

(2) In particular, for any z G C \ Vs^r{Pri') we have 

\\(z - Pn^y^Wsv < iiii:^r-"i + — < 1036.693385. 
1 — r 2ei 

Recall that ei = — ''"^^"f — . 

(3) Now we go back to the problem of finding the size of a hole which guarantees 
the existence of a Tff-accim with the desired escape rate. Here a — 3ao 
and all we have to do is to feed the estimate on Hg.riPri) obtained in Step 
2, together with the new 712, in the formula of Sq to obtain that 

i^ryha > 0.00001216687545. 

Hence, we can deduce that mesh(77') — 10^'"' will do the job; i.e., X{H) G 
(0,-^ x 10^^] =^ Th admits an accim /i with escape rate — Ine/f < 
-ln(39/40). 
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Upper bound on Hs,r 


63.73181657 


1036.693385 


ni 
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40/39 


40/39 


n2 
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Lower bound on (2r)~^eo 


0.0001763820641 


0.00001216687545 


Loop I 


Fail: reduce e 


Pass 


Loop II 




Pass 


Output I 




Ecom = 10'\(5com = 1/41 


Output II 




MH) G (0, ^ X 10"^] =^ Th admits an accim ^l 
with escape rate — Inen < — ln(39/40) 



Table 2. The output of Algorithm O for / = 1/40 



7. The effect of the position of a hole 

The results of the previous section give upper bounds on the escape rate that are 
uniform for a given size of hole, independent of the position of the hole. However, it 
has been observed already in [4J that the position of the hole can affect the escape 
rate; i.e., given a map T and two holes Hi, with X{Hi) — X(H2), it may happen 
that the escape through Hi, say, may be bigger than the escape rate through H2- 
For example, define the sets 

Per(i7,) ^{p: peN s.t. for some x G Hi, TP{x) = x, and TP-^{x) ^ x}; i ^ 1, 2, 
For certain maps, if 

Minimumjp G Per(_ffi)} < Minimum {p £ Per(_ff2)} 
then the escape rate through Hi will be smaller than the escape rate through H2 ■ 

In [llj Keller and Liverani obtained precise asymptotic information about the 
effect of the location of the hole. Roughly speaking, for a system of holes shrinking 
to a single point, the rate of decay of escape rate depends on two things: the value 
of the invariant density of the map T at the point the holes shrink to, and whether 
or not this point is periodic. We now state a version of this result and will discuss in 
the next subsection, in a smooth setting, how a combination of Algorithm lS.ll with 
the proper modification of P^, can be used with this theorem when the formula 
of the invariant density of T cannot be found explicitly. When smoothness is not 
assumed, as in this paper, obtaining asymptotics for the escape rate relative to the 
size of the hole appears to be an open problem. 

Theorem 7.1. [llj Let T be piecewise on a finite partition of [0, 1] and assume 
it is piecewise expanding and mixing. Let {iJ^} be a sequence of holes such that 
Hf^ D H^i , for < k' < K, with Hq — {y} for some point y £ [0, 1] which is a 
point of continuity of both T and f* , f* is the invariant density of T . Let Th^ be 
a perturbation of T into a map with a hole. Assume that inf f*j^ > 0. For X{Hi^) 

sufficiently smal^^ we have: 



can of course quantify what we mean by sufficiently small using Algorithm l5.ll 
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(1) If y is non-periodic then liiUK^o = f*{v)- 

(2) Ify is periodic with period p then lim^^o ^\(h") = f*{v) (l - \{Tpy(y)\ ) ■ 

7.1. circle maps. Theorem 17.11 requires the knowledge of the value of the 
invariant density /*; in particular, its value at the point y. Unfortunately, the 
approximate invariant density which is obtained by Ulam's method in Algorithm 
IS.ll does not provide a pointwise approximation of /*. However, in a smooth setting, 
one can modify Ulam's scheme, and the function spaces where P and act, to 
obtain rigorous approximation of ||/* — /,,||oo, where denotes the fixed function 
for the modified finite rank operator P^. Of course our main theoretical tool, 
Theorem 12.21 will need to be modified. In [lOj an abstract version of Theorem 12.21 
is proved. The result requires bounded linear operators Pi and P2 acting on two 
abstract Banach spaces whose norms || • || and | • | satisfy | • | < |1 • ||, and the unit 
ball of II • II is I • |-compact. Thus, Theorem 12.21 is a particular application of the 
general result of [TU]. When T is a circle map, it is well known that its invariant 
density /* is C^. Thus, instead of and BV , one can study the action of P and 
a smooth versio of on the spaces W^'^ and W^'^ with norms 

ll/lk^.^-||/lli + ll./'lti 

ll/IU.,. = ||/||i + |l/'||i + ||/"||i 

respectively. A common Lasota-Yorke inequality of P and Pri in this setting is given 
by: for / e W^^'^ and n e N we have 

||P,/||;^.. =a2"||/||^^i,.+5||/||H.i.i, 
where _B > which depends on T only. Using this setting, one obtains the estimate 

II/* -/»?l|wi.i < C'-mesh(?7). 

For more details and for a proof of the above Lasota-Yorke inequality we refer to 
Section 10.2 of [ig. 

In a setting like this, one can then repeat Algorithm 15.11 with the smooth version 
of Pri and obtain the following reformulation of Theorem [73] to a setting where the 
invariant density /* is a priori unknown. Note that for / e W^'"^, \\f\\oo < ||/||ivi i- 

Theorem 7.2. Let T be a C'^ circle map. Let {H,^} be a sequence of holes such 
that iJ^ ^ Hk', for < k' < k, with Hq = {y} for some point y £ [0, 1]. Let f* 
be the invariant density of T , and Th^ be a perturbation of T into a map with a 
hol^^. Let £ = mesh(rj). 3 a constant C = C{Pr^) such that for \{Hk) G {Q,Tecom\, 
we have: 

(1) Ify is non-periodic then 

fviy) -C-e< lini < f^{y) + C-e. 



^^For instance one can use a piecewise linear approximation method |6]. 

""^^inf^gp^l] /* > for circle maps. See [12] or [IT]. Thus, the assumption inf /p^ > is 
automatically satisfied for such maps. 
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(2) // y is periodic with period p then 

(«=')-C-.)(i-^)<l»i^< 

Proof. We only give a sketch of the proof. Suppose that we have used Algorithm 
15. II with the proper modification of Pri and the function spaces. Then the invariant 
density, which is a byproduct of the algorithm, would provide the following estimate: 

ii/„-riioo <c.£. 

Consequently, for any y G [0, 1], we have 

(7.1) \fn{y)-nv)\<C-e. 

□ 

Thus, the proof follows by using (|7.ip and Theorem 17. II 

Remark 7.3. All the constants which are hiding in the computation of C = C^Prj) 
can be rigorously computed using Theorem 12.21 with the spaces W^'^ and W^'^. It 
should be pointed out that these constants cannot be computed if one attempts to 
do this a approximation in the L^, BV framework. This is because the estimates 
will depend on (/*)" which is a priori unknown. 

Remark 7.4. For Lasota-Yorke maps, a result similar to Theorem 17.21 is not 
obvious at all. The problem for Lasota-Yorke maps involves two issues: 

(1) The invariant density /* is . This means that the framework of VF^^^, 
W^ "^ cannot be used. 

(2) If one uses the function spaces BV and L^, then to the best of our knowl- 
edge, only the original Ulam method will fit in this setting. The problem 
with Ulam's method is that it provides only good estimates in the norm 
II/* - /r,||i = C • elnl/e. However, typically, ||/* - /^Wbv 0. 

Our last comment on this is that providing a scheme for Lasota-Yorke maps to 
obtain a result similar to that of Theorem 17.21 would be an interesting problem. 
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